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Using the k-p theory and first-principles simulations, we report that applying a moderate pressure 
(> 0.6 GPa) on black phosphorus can diminish its band gap and produce one-dimensional and 
even two-dimensional (2D) Dirac cones, distinguishing this material for use in novel non-compound 
topological insulators. Similar to topological insulators, these 2D Dirac cones result from two 
competing mechanisms: the unique linear band dispersion tends to open a gap via a “pseudo spin- 
orbit” coupling, while the band symmetry requirements preserve the material’s gapless spectrum. 
Moreover, these unique Dirac cones are bulk states that do not require time-reversal symmetry, thus 
they are robust even in the presence of surface or magnetic perturbations. Ultimately, we show that 
our predictions can be detected by the material’s unusual Landau levels. 

PACS numbers: 


To date Dirac materials, such as graphene [1 i, 
topological insulators (TIs) or spin Hall insulators 
(sms) @, 0, topological crystalline insulators (TCIs) 
d, 0 and Weyl semimetals (WSs) [13, [III , have drawn 
tremendous research interest. Graphene hosts two- 
dimensional (2D) Dirac fermions [H, HI]; TIs are mate¬ 
rials with a bulk energy gap but possess gapless 2D (ID) 
Dirac surface states that are protected by time-reversal 
symmetry [13, [i3| ! TCIs exhibit metallic surface states 
protected by the mirror symmetry of the crystal 0; WSs, 
which are protected by crystal symmetries, exhibit en¬ 
ergy overlaps that occur only at a set of isolated points 
in momentum space where the linearly dispersed bands 
cross at the Fermi level 0- The impact and implications 
of these materials are difficult to overstate. The discov¬ 
ery of the anomalous Hall effect fla UTII . the fractional 
Hall effec t 113. chiral anomaly [3, [^, and Majorana 
fermions |2lL [22| , to name a few, have inspired enormous 
efforts in search for new Dirac materials. 

Recently a new type of 2D semiconductor, few-layer 
black phosph orus (BP) (phosphorene), has been fabri¬ 
cated |23l425j . It exhibits promising carrier mobilities 
and unique anisotropic transport, optical, and thermal 
properties (23 - E3 |. In particular, few-layer BP posses a 
unique band structure, whose dispersion is nearly linear 
along the armchair direction but parabolic along other 
directions [13, [Hi. An obvious idea follows from this 
unique structure: if one can significantly reduce the band 
gap or even achieve the band inversion, novel features 
may result, e.g., the formation of Dirac cones. However, 
the band gap of monolayer and few-layer BP is signifi¬ 
cant (up to 2.2 eV) [13, [ 23 , mi , making it impractical to 
close. Alternatively, bulk black phosphorus has a much 
smaller band gap, only around 300 meV [ 13 , [ 131 , which 
is more tractable. 

In this work, we show that the k ■ p Hamiltonian of 
bulk BP can be written as the free-electron-gas Hamil¬ 
tonian with off-diagonal terms for describing interband 
interactions, which are of the same form as the Rasshba 
and Dressehaus spin-orbit coupling (SOC). Thus we call 
this off-diagonal interaction the ” pseudo spin-orbit” cou- 
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FIG. 1: (Golor online) (a) Top view of the atomic structure 
of bulk BP. The unit cell and lattice vectors are marked. The 
zigzag and armchair directions are illustrated as well, (b) Side 
view of BP. (c)The first BZ. 


pling (PSOC). If the band gap is diminished and the band 
inversion is realized, PSOC tends to open a gap in the 
energy spectrum. However, similar to TIs, the band sym¬ 
metry requires band connections. Competition between 
PSOC and band symmetry causes the Fermi surface to 
undergo a topological transition from a sphere to a ring 
of 2D Dirac cones. Using first-principles simulations, we 
show that pressure can be an efficient tool to produce 
these 2D Dirac cones and even ID Dirac cones, while 
also tuning the characteristic Fermi velocity. Interest¬ 
ingly, because our predicted Dirac cone feature in com¬ 
pressed BP does not require the time reversal symmetry, 
it is not destroyed by magnetic fields. Thus we provide 
a means for detecting this unique electronic structure by 
demonstrating the unusual energy scaling law exhibited 
by Landau levels. 

The atomic structure of bulk BP is plotted in Figures 
1 (a) and (b). Here we use a simple orthorhombic con¬ 
ventional unit cell, which is twice the size of the primi¬ 
tive cell. The corresponding first Brillouin Zone (BZ) is 
shown in Figure 1 (c). The double-cell procedure sim- 
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FIG. 2: (Color online)(a) and (b) are the band structures of 
intrinsic BP. The red dots are from the HSE06 simulation and 
the blue line is from the k -p calculations, (c) and (d) are the 
band structures calculated from the k ■ p calculations under 
the band-inversion condition. The Fermi level is set to be 
zero. 

ply folds the first BZ and does not affect any essential 
physics; this is convenient for presenting the topology 
of the Dirac cones. The bulk BP atomic structures is 
fully relaxed according to the forces and stresses calcu¬ 
lated using density functional theory (DFT) with the 
Perdew, Burke and Ernzerhof (PBE) functional associ¬ 
ated with van der Waals corrections [13, Because 
DFT is known to underestimate the band gap , we use 
a hybrid functional theory (HFT), HSE06 [3a js^l, to 
calculate electronic band structures, which is in good 
agreement with the GW calculations and experimental 
measurements [13, dll of bulk BP. 

We begin with the widely used k ■ p model to capture 
the essential BP’s electronic structure. The D 2 h point 
group invariance allows us to describe BP using a two- 
band model similar to the one used for monolayer BP 
[13, 113]. Thus the low-energy effective Hamiltonian is 
given by 

f Ec + akl+ Pkl + 'ykl Vfky \ 

Vfky Ey - Xkl - pkl - vkl, ) 

( 1 ) 

where kx^ ky, and kz are components of the crystal mo¬ 
mentum along the x, y, and z directions, respectively. 
The conduction band and the valance band energies, 
Ec = 0.15 eV and Ey = —0.15 eV, are set to fit the 
bulk band gap (300 meV). a, ,5, 7 , A, and v are fitted 



FIG. 3: (Golor online) The 3D band structures of BP of the 
X-Y plane (a), the Y-Z plane (b), and the X-Z plane (c). The 
band overlap in (c) forms a ring structure, (d) The evolution 
of the band structure from the X-Z plane to the X-Y plane 
via changing of the tilted angle. 


parameters that describe the band edge curvature, and 
Vf is the Fermi velocity of the nearly linear band disper¬ 
sion along the armchair (y) direction. By appropriately 
choosing these parameters, as shown in Figures 2 (a) and 
(b), the k ■ p result nearly perfectly matches that of the 
first-principles HSE06 calculations of intrinsic BP. 

The special characteristics of Eq. 1 are the off-diagonal 
terms of the effective Hamiltonian, which produce the 
unique linear dispersion along the armchair (y) direction, 
as shown in Figure 2 (a). In fact. This off-diagonal inter¬ 
action can be split into two terms 

^off = Husoc + Husoc 

= 1/2 Vf{k^<jy + kya^) - 1/2 Vf{kz,ay - kya^), 

where cr^ and ay are Pauli matrices. Interestingly, these 
two terms {Hdsoc and ffff.s’or?) have the forms of the 
Dressehaus and Rashba SOC [40], respectively. Realis¬ 
tically, Dressehaus and Rashba SOC cannot exist in in¬ 
trinsic bulk BP, because its lattice structure possesses 
the inversion symmetry. Hence we call these interactions 
PSOC. Because of their similar formalisms, we expect 
equivalent ”SOC effects” to appear in bulk BP. In par¬ 
ticular, considering the dramatic effect of SOC in TIs, 
it is natural to expect that diminishing or inverting the 
band gap in bulk BP will produce novel electronic char¬ 
acteristics. 

It is easy to mimic the band inversion in the k ■ p 
model; one simply assigns the valence band edge energy 
Ey a higher value than the conduction band edge energy 
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FIG. 4: (Color online)(a) The HSE06-calculated band struc¬ 
ture of BP under the critical pressure (0.6 GPa). (b) The 
HSE06-calculated band structure of BP under the 0.9 GPa 
pressure, (c) The Anisotropic fermi velocity of electrons and 
holes of 2D Dirac cones varied with the pressure. 


Ec in Eq. 1. Surprisingly, instead of producing simple 
band overlaps, we observe an unusual topological transi¬ 
tion at the Fermi surface. In Figures 2 (c) and (d), the 
band crossing at the F point is opened by the PSOC, 
while the band crossing is preserved at the K point in 
the r — X direction. As a result, 2D graphene-like Dirac 
cones are formed and the Schematic 3D plots of the elec¬ 
tronic structure are presented in Figures 3 (a), (b), and 
(c) for three specific cutting planes, X-Y, Y-Z and X-Z, 
respectively. The linear band dispersions and Dirac cones 
are observed in both X-Y and Y-Z planes and a ring-like 
band crossing occurs in the X-Z plane. Actually, we can 
regard the Dirac cones observed in Figure 3 (a) and (b) 
are those points marked in the ring structures of Figure 
3 (c). In this sense, any point on the ring of Figure 3 
(c) corresponds to a Dirac cone in other 2D planes. In 
Figure 3 (d), we show that this ring-like band crossing 
continuously evolves to a pair of Dirac cones by varying 
the tilting angles from the X-Z plane to the Y-Z plane. 
In other words, the set of tilting angles that do not ex¬ 
hibit these Dirac cones is of measure zero, and thus our 
findings are robust. 

The mechanisms responsible for forming the ring of 
Dirac cones are similar to those in TIs, but there are some 
crucial differences. First, similar to the band inversion 
and gap-opening that occurs in TIs due to SOC interac¬ 


tions, the gap-opening at the T point of Figure 2 (c) from 
the PSOC of states residing on the linear band along the 
r — Y direction. Second, the gap of BP is not completely 
opened; the valence and conduction bands still touch at 
the K point (Figure 2 (c)) or by the ring (Figure 3 (c)). 
This gapless feature is guaranteed by a symmetry con¬ 
straint based on two facts: the bottom conduction band 
is formed by what was originally top valence band prior 
to the band inversion; the parity symmetries of the top 
of the valence band and the bottom of the conduction 
band are inverse to each other, which is evidenced by the 
bright dipole optical transitions [2^. Therefore, as shown 
in Figures 2 (c) and 3 (c), contact between the conduc¬ 
tion and valence band is required in order to provide a 
continuous symmetry evolution. This mechanism is sim¬ 
ilar to the one responsible for forming gapless surface 
states in TIs [4l|. However, there are no surface states 
in our calculated bulk BP, thus it is the bulk valence and 
conduction bands that must touch. As a result, these 
two factors force the formation Dirac cones. Moreover, 
as proposed by a recent work [42j, further breaking the 
inversion symmetry and including the spin-orbit coupling 
in few-layer BP will open the bulk band gap and produce 
surface topological states. Finally, another special char¬ 
acteristic of these bulk Dirac cone states in BP is that 
they are robust even in the face of broken time reversal 
symmetry. Therefore, these Dirac cones can uniquely be 
observed in experiments that employ magnetic fields and 
are potentially useful for magnetic devices. 

Following the above analysis of the k ■ p model and be¬ 
yond 2D Dirac cones, 3D Weyl points could be achieved 
if an extra direction, such as the x or z, also exhibits the 
off-diagonal linear interband interaction (PSOC) term. 
Therefore, one may search for WSs by looking for re¬ 
alizing the band inversion of the materials that exhibit 
nearly linear dispersions along two orthogonal directions. 

Beyond the k ■ p model, an obvious question is how 
one can practically reduce the band gap or achieve band 
inversion in BP. Towards this end, we turn to the first- 
principles simulations for quantitative answers. Our 
HSE06 calculations show that applying an external pres¬ 
sure along the y (armchair) or x (zigzag) direction pro¬ 
duces promising effects. Note that, BP is much softer 
along the y (armchair) direction because of its anisotropic 
atomic structure [i^. We find that applying a pressure 
of 0.75 GPa can shrink the lattice constant along the y 
direction by 4% but only by 1% and 2% for the x and z 
directions, respectively. Therefore, we will consider the 
changes made by applying pressure along the y direc¬ 
tion because it can change the atomic structure and cor¬ 
responding electronic structure more significantly via a 
relatively small pressure. 

First, as shown in Figure 4 (a), an applied critical pres¬ 
sure (0.6 GPa) diminishes the band gap and a ID Dirac 
cone is produced along the P — A direction. This pres¬ 
sure is well within the capability of current high-pressure 
experiments [3, . 

Second, larger pressures create a band inversion. Ex- 
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FIG. 5: (Color online) (a) The energy spectrum of Landau 
levels (index n) of intrinsic BP with magnetic field; (b) that 
under the critical pressure (0.6 GPa), in which the ID Dirac 
cone is formed; (c) That under the pressure of 0.9 GPa, in 
which 2D Dirac cones are formed. 

citingly, similar to the k -p theory predictions, a pressure 
of 0.9 GPa yields a band crossing at the P point and a 
pair of 2D Dirac cones are formed, as shown in Figure 4 
(b). The corresponding density of states (DOS) is shown 
in Figure 4 (c). We can clearly see that a linear DOS 
in the low-energy regime, which is similar to graphene. 
From these first-principles results, we also find that that 
the Fermi velocities of compressed BP is smaller than 
that of graphene and that these generated Dirac cones 
are highly anisotropic in momentum space. Obviously, 
the Fermi velocities must be anisotropic as well and can 
be efficiently tuned by varying the applied pressure, as 
concluded in Figure 4 (d). 

We propose that the evolution of the electronic struc¬ 
ture from conventional semiconductors to one possessing 
ID Dirac cones and 2D Dirac cones can be observed by 
measuring the energy spectrum of Landau levels; this 
methodolo^ is widely used for identifying Dirac band 
structures Q- For example, considering the X-Y plane, 
before reaching the critical point for closing the band 
gap, the energy spectrum of the Landau levels is that of 
a typical semiconductor, i.e., it is linear with the level 
index, as shown in Figure 5 (a). Here the spectrum is 
given by 

E = Ec -l- (n -|- l/2)?iWc 

E = Ey — (n + l/2)hujy, 

where the cyclotron frequencies of electrons and holes are 
Wc = eB/y/rricxmcy and tUy = eB/^my^niyy. rricx, mcy, 
myx and ruyy are effective mass of electrons and holes 
along the X and Y directions, respectively. Importantly, 


there is no zero-energy solution for this case. 

At the critical pressure, a perfectly ID linear band dis¬ 
persion along the F — Y direction simplifies Eq. 1 by set¬ 
ting the parabolic term /3 and /i to be zero. The effective 
Hamiltonian with the magnetic field introduced by the 
Landau gauge ^ = {—By, 0, 0) is 

-^^{hkx - eByf(f)A - ihvf-^(j)B = E(j)A, 

Zrrir dif , 

1 d 

- - {hkx - eBy)'^(j)B - ihvf—(j)A = E(j)B, 

Zrriy ay 

where vj = 3.5 x 10®to/s is the Fermi velocity of the ID 
dirac cone obtained from HSE06 simulations. The corre¬ 
sponding energy spectrum of Landau Levels in magnetic 
field is presented in Figure 5 (b), which is deviated from 
the linear relation of Figure 5 (a). Finally, for higher 
pressures that form 2D Dirac cones, the Landau levels 
follows the well-known square-root dependence [1, [l^ . 
En oc y/n , and exhibits a unique zero-energy solution, 
as shown in Figure 5 (c). 

In conclusion, we predict that applying a moderate 
pressure (> 0.6 GPa) on bulk BP can dramatically re¬ 
shape its band structure. Its band gape will close and 
band inversion will induce a topological transition of the 
Fermi surface via the so-called PSOC interaction. Un¬ 
der different pressures, the band dispersions of bulk BP 
will evolve from that of a semiconductor with a parabolic 
band dispersions to a semimetal with ID or 2D Dirac 
cones that exhibit tunable anisotropic Fermi velocities. 
Gompressed BP thus may be a TI that uniquely inte¬ 
grates both massless Dirac Fermions and massive clas¬ 
sical carriers within one material. Furthermore, these 
bulk Dirac-cone states and robust even under broken 
time reversal symmetry. Therefore, these unique elec¬ 
tronic structures can be identified by Landau levels and 
we have provided the corresponding energy spectra. 
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